Nuclear surface studies with antiprotonic atom X-rays 
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The recent and older level shifts and widths in p atoms are analyzed. The results are fitted by an 
antiproton-nucleus optical potential with two basic complex strength parameters. These parameters 
are related to average S and P wave scattering parameters in the subthreshold energy region. A 
fair consistency of the X-ray data for all Z values, stopping p data and the A*' A*' scattering data has 
been achieved. The determination of neutron density profiles at the nuclear surface is undertaken, 
and the determination of the neutron Rrms radii is attempted. Uncertainties due to the input data 
and the procedure are discussed. 
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I. INTRODUCTION 



Recent CERN/LEAR measurements of X-ray transitions in antiprotonic atoms provide new data on the level shifts 
and widths in heavy systems These results complement older CERN, BNL and other [3, [1,0] studies in medium 
Z atoms and fairly recent experiments with hydrogen, deuterium and helium 0, SSQ- The main distinction of the 
^ ' new results from those obtained earlier in light nuclei is a much broader data basis, extended to the region of neutron 
^\ , excess and deformed nuclei. With the enlarged data one can attenipt an extraction of a phenomenological p optical 
OA potential that contains separate strengths of and np interactions It allows to study properties of the neutron 

density distributions at the nuclear surface. As stressed on many occasions the p atoms are well suited for such studies. 
A problem is met on this way as the best description of the atomic data indicate the pp and rip absorption rates to be 
roughly the same. This result is in conflict with the chamber studies of low energy p annihilation [III, [I^l- The latter 
, experiments indicate the pp absorption rate to be twice as fast as the rip one. This discrepancy has to be understood 
' and we indicate the solution, 
"•pi ', The content of this paper is as follows. In section II, the optical potential is constructed for p bound into atomic 
' I ' states. It consists of two terms, the central and the gradient one. The strengths of these terms are phenomenological 
O [ parameters which describe the S and P wave NN scattering amplitudes extrapolated below the NN threshold. The 
5 ■ easiest way to learn about such amplitudes are the simple systems ^'^H— p and '^'^He— p. Here, we rely on the absorptive 
amplitudes extracted from these atoms. Additional constraints follow from the chamber experiments [ill, [l^. Some 
degree of unification of these data is obtained in terms of NN interaction potentials, in particular the recent updates 
of the Paris potential [l3| . This allows for an approximate separation of the S and P wave absorption which is 
^ , implemented into the optical potential. 

• The X-ray data from the four lightest atoms indicate the existence of two quasi-bound NN states. One state occurs 
in an 5 wave. It finds strong support in the J/tp — » pp ^ decays which allow to pinpoint its quantum numbers, [l^ . 
The second quasi-bound state occurs in a P wave, [l5|, and so far has no independent confirmation. Both states are 
reproduced by the Paris model, [l^, and this facilitates our discussion of the p atoms. 

Second question studied in Section II is the sensitivity of atomic levels and level widths to the nuclear densities. 
A significant dependence on the input charge densities is found. We attempt a model independent parametrization 
of the neutron densities and argue that p atoms test nuclear surfaces but not the single particle asymptotic density 
regions. 

The optical potential parameters are found via the best fit procedure to 117 atomic X-ray data. The fit is improved 
considerably by the effect of NN S wave quasi-bound state. The other, P wave, state is of no importance but it can 
explain anomalies observed in p capture on loosely bound nucleons [16,] . 

In a number of cases the atomic level widths may be rather precisely measured for two orbits per atom, the "lower" 
and the "upper" one. Such widths are useful to study properties of the nuclear surface. Atomic level shifts are 
less accurate and more difficult to understand. These difhculties reflect to a large extent the uncertainties in the 
understanding of basic NN interactions. 

In section III, the Rrms radii of the neutron density distributions - Rrms - are extracted for several isotopes of 
Ca, Zr, Sn, Te and Pb. This is done on the basis of X-ray data and radiochemical measurements [13, [iM [isl - [20II 
which test the ratios of neutron and proton densities in the region even more peripheral than that for the X-rays [211 ] . 
The radiochemical data determine the rate of pn capture relative to pp capture at very large distances. These are 
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discussed as a separate issue, since the nuclear region involved in this process is located about 1 fm farther away from 
the region tested by the atomic X-rays. Most of the effort is devoted to the evaluation of uncertainties involved in 
this method of density determination. 

The appendix discusses some details of the gradient potential. 

The present publication addresses also three more specific questions: 

• The optical potential for antiprotons involves uncertain quantities: the range and strengths of the NN interaction 
and nuclear densities at large distances. Errors due to these uncertainties on the atomic level widths are 
evaluated. 

• The experimental level widths in heavy atoms are determined by high moments of the neutron density. However, 
the nuclear structure physics is more interested in the low moments, in particular in the rms radius. How well 
could we determine the latter? 

• Can antiprotonic data distinguish a neutron skin from a neutron halo as defined in Rcf. '22\1. 

A number of phenomcnological optical potentials have been fitted to the p X-ray data. Recent results may be found 
m Refs. ^E^. The present work is different in several aspects: 

• We include recoil effects in the P wave NN interactions. 

• The constrains from Z = 1, Z = 2 atoms, p absorption in flight and NN potential models are accounted for. 

• The data set is larger. 

The constraints allow to obtain the absorptive optical potential parameters close to those expected from the ^A'^ 
scattering data. In consequence a good fit to the atomic X-ray data is obtained. 

II. THE OPTICAL POTENTIAL 

Atomic energy levels of antiprotons are determined essentially by the Coulomb and fine structure interactions. In 
addition, the p- nucleus interactions generate level shifts e and p annihilation generates level broadenings T. Both 
effects may be conveniently described by a complex nuclear optical potential. The standard potentials, well tested for 
TT atoms [!|| , are composed of two terms 

t/°P* = Ep,„ [Vs{r)+WVp{r)V]^Vs + Vp, (1) 

with the sum extending over protons and neutrons. Both terms, the local and the gradient Vp are expected to 
have a folded form 

Vs,p{r) = as.p / du gs.p{u)p{r - u), (2) 

where ^iffj^ is the NN reduced mass and p is the nuclear density. Here as is expected to resemble the spin averaged S 
wave scattering length and ap the P wave scattering volume. Two profile functions gs,p, normalized by / dugs,p{u) = 
1, reflect the range of interactions. The NN annihilation radius is expected in the range 0.8 — 1.0 fm. Such values 
follow from phenomenological or quark models of A^A^ interactions. The effective annihilation radius in models with 
much shorter annihilation potentials is similar [l^]. On the other hand, the ranges involved in Re Vs and Re Vp 
may be different. In this work, the same interaction range is assumed for all the components, and it is left as a free 
parameter. A Gaussian profile g is used. 

The form of given in eq. ([T]) is related to the parametrization of the low energy scattering amplitudes 

/ = as + Spp'ap (3) 

where p, p' are the relative momenta of the colliding particles before and after the collision. An important distinction 
between NN and ttA'' cases is that in the A^A^ collisions the nucleon recoil effect is important and the relative 
momentum is p = (pjv — Pjv)/2- It involves the nucleon pat and antiproton p^ momenta in equal proportions. As a 
consequence, gradients over nucleon wave functions arise in the optical potential. One needs a formula that generalizes 
eq.([T|). The scattering amplitude ([3]) leads now to a new form of the gradient potential 

Voir) = ^ap^I]„ / drV:(r')(^iv - ^jv)/p(r " r'){^ n " ^Jv)^c(r') (4) 
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where </3q, are the nucleon wave functions. The summation over nucleon states a and the integration over nuclear 
coordinates r' are to be performed. This leads to a three-term expression for Vq 

Vg = Vp + Vn + V™„ (5) 

and 

l/°P* = Sp,„ [Vs + Vg]- (6) 

The first term in Eq. ([5|) corresponds to the standard gradient potential, Vn is due to the gradients over nucleon 
functions, while Vmix follows from mixed nucleon and antiproton derivatives. In the states studied in experiments 
the nucleon dependent part Vn + Vmix contributes about half of the Vg strength, and amounts to a quarter of the 
total yp*. On a phenomenological level these two terms could be included into the Vg potential term. However, 
the difference arises when one attempts to relate as and ap to the scattering data. A special effect comes from the 
dependence of these terms on the state of the nucleus, in particular on the angular momentum in the valence shells. 
In addition the mixing term may induce some nucleon-antinucleon correlations in odd-A nuclei. A more detailed 
discussion of the gradient terms is given in Appendix A, where some approximations are also introduced. 

A. The choice of nuclear charge densities 

The nuclear charge densities are well determined in the region of the nuclear surface between c — 2a and c + 2a, 
where c is the half density radius and a is the surface diffuseness. These are the results of electron scattering and 
muonic atom experiments psl |23 |. However, the antiprotonic atoms involve also lower nuclear densities. In Fig.l 
are shown the average radii of p absorption in the "lower" atomic orbits. At those radii the nuclear charge densities 
amount to 5% of the central density. For the "upper" levels these radii are larger by 0.2 — 0.4 fm and involve charge 
densities smaller by a factor of 2. The absorption regions are localized in nuclear layers of about 3 fm radial depth. 
In these regions, the charge densities are not well known. To indicate the uncertainty, let us compare the relevant 
moments of several charge distributions. Atomic level widths in very high I orbits are given by expectation values 
r/2 ~< n,l\V°P*\n,l > where / is the angular momentum and n the principal quantum number of an atomic state. 
For large Bohr radii and weak nuclear absorption, the widths are proportional to < r^' >, the 21— th moment of the 
nuclear density distribution. This reflects the dominant effect of high centrifugal barriers. In the states actually tested 
the finer details of the p wave functions are significant and the dominant moments are < r^'~^ > and < r^'~^ >. 
Several relevant moments of nuclear charge density are compared in Table IH For a given isotope, all charge density 
profiles yield essentially the same rms values. However, the moments of interest may differ by up to 30% in the highest 
observable I states. In particular the many-parameter multi-Gaussian density parametrization offer much shorter tails. 
Other extreme cases are given by some three-parameter Fermi distributions which often generate unrealistically low 
(and even negative) densities at very large distances. To select the "best" charge profile we follow the "averaging 
procedure" outlined in Ref. [2^ for hyperonic atoms. Thus, the lowest moments < >, < r'' >, < > are compared 
for several available charge profiles (Fermi, Gaussian, multi-Gaussian) and an average density in the sense of average 
moments is extracted. In the cases of Al, S, Ca, Pb, studied in Ref. [1^ one always finds the average to be the 
closest to the profiles given in the most recent compilation by Fricke et al The same is found in all Ni, Zr and 

Sn isotopes studied here. We use the parameters of Ref. [24] as the basis of our calculations. The only exception is 
^^O, where the sum of Gaussians and the deformed nuclei U, Th, where monopole density component were used [231]. 
Later, in specific cases, the comparison with other densities is presented. 

B. The parametrization of nuclear densities 

To understand the p atomic data, one needs a reliable extrapolation of the densities to very large distances. The 
related problems are visualized in Fig. 2. For two Ca isotopes the densities were calculated with a HFB-SkP model [26| , 
and fitted at large distances by a two-parameter Fermi profile 

p(r) = l/(l + exp[(r-c)/a(r)]. (7) 

Instead of a constant diffuseness one has to introduce certain functions a(r), which are plotted in Fig. 2. The 
dependence on radial distance is rather distinct. There are clear advantages of the 2pF profile, but with a constant a 
it does not reproduce the density in far away regions. That result is also fairly model independent and the shapes of 
a{r) in Fig. 2 indicate a certain degree of universality. At distances c + 3a<r<c + 6a one finds essentially the same 
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FIG. 1: Dots - the average antiproton absorption radii < Rl > (weighted by t/jatom I™ V°'''r^ ) calculated for the "lower" 
atomic level, left scale in fm. Squares - the charge densities at these points, given on the right scale as % of the central charge 
densities. The bottom scale - atomic numbers Z. Calculations are done with densities from Ref. i24i1. 



slope of a(r) for a wide range of nuclei in Hartree-Fock, Hartree-Fock-Bogolubov [2l|, [2y| and Relativistic Mean Field 
calculations [2?|. However, at radii beyond this region, different nuclear models may produce different behavior. 

The average annihilation radii ra in Ca are marked by arrows. In the lower and upper states one has ra = 4.9 fm 
and 5.3 fm respectively and a(r) is seen to be fairly stable at these ranges. However, the radiochemical data involve 
larger radii from « 5.1 fm up to r « 7.5 fm ^2ll|. The diffuseness parameter a(r) is seen to fall down in this region. 
The difficulty involved in the radial dependence of a{r) is moderated by the optical potential which involves folded 
densities. With the folding range of « 1 fm the corresponding downfall of ajoidedif) is pushed away by « 0.5 fm. 

In order to understand a{r) the first question to answer is where the single particle asymptotic limit is reached and 
what is its form. For a single neutron of angular momentum L bound in an external potential the asymptotic wave 
function is described by 



exp(— rfca) 



(8) 



where = \/2mpfEa, Ea is the neutron binding energy and Wl is a polynomial given by the Hankel function of 
order L (see e.g. ^]). This wave function allows to scale the asymptotic density p{r) — (pa{r)'^. The ratio of two 
densities p{r)/ p{ro), taken at two points and r, allows to find 
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1 



a(r) a(ro) 



WL[r)ro r 



(9) 



This formula reproduces the shapes of a(r) in the region of r « 5 — 8 fm with ~ 1.25 fm~^ which corresponds to 
Ea ~ 30MeV. On the other hand, the separation energies are respectively 16 MeV (^°Ca) and 10 MeV (^^Ca). Thus 
the "true" single particle asymptotic density given by valence neutrons begins farther away. The second observation 
is that the effective fc^, which reproduces Fig. 2, corresponds to bindings larger than the average binding weighted by 
contributions of neutron orbitals to the total neutron density. The behavior of a nucleon at large distances involves 
virtual excitations of the residual system and a sizable fraction of ka is due to the nuclear correlations. Contrary to 
a frequently expressed belief, the radii tested by antiprotons do not represent the far away distances of single particle 
asymptotic wave functions, and in this sense are more interesting for the nuclear structure research. 

In practical terms the shortcomings of constant a can be corrected by the a{r) given by eq. ^ with ka ~ 1.25 
fm^^. Because of the large value of ka the result is in practice independent of the value of L. Relation ^ is used 
here to interpolate between two values of the diffuseness parameter. The initial value a(ro) taken at = c + 3a is due 
to the full complexity of the nuclear structure which includes the average field and effects of nuclear correlations. It 
is kept as a free parameter, the a, to be determined from experiments and reproduced by models. Such a procedure 
serves only as a guiding principle for the best correlation of the atomic X-ray data and the radiochemical data. With 
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TABLE L Moments < r^' >= J drp(r)r^7 J drp(r) of the charge density profiles in 10^'"^ fm^' units. The densities based 
on jj, atom data -first column, electron scattering data - other columns: 2pF - two parameter Fermi, HO - harmonic oscillator, 
SOG - multigaussian. 



Atom 


moment 


2pF(^) [24] 


2pF(e) [23] 


HO(e) [23] 


SOG(e) [23] 


16(3 


rms 


2.72 


2.72 


2.72 


2.71 




< > 


9.0 


8.7 


8.2 


8.4 




< r'^ > 


1.6 


1.5 


1.2 


1.3 


32s 


rms 


3.28 


3.25 




3.26 




< r* > 


17.2 


16.9 




1.69 




< r« > 


3.75 


3.81 




3.82 




rms 


3.79 


3.78 




3.77 




< r« > 


7.5 


7.6 




6.7 




<r^ > 


2.6 


2.5 




1.9 


"'^Sn 


rms 


4.64 


4.65 




4.63 




< r** > 


8.6 


9.8 




7.9 




< r"' > 


4.2 


5.1 




3.7 




<r" > 


2.5 


3.3 




1.7 


i^^Sn 


rms 


4.68 


4.66 




4.69 




< > 


9.3 


9.5 




8.4 




< r"' > 


4.6 


4.8 




3.7 




<r'^ > 


2.6 


2.9 




1.8 


208 


rms 


5.504 


5.520 




5.503 




<r'^ > 


11 


13 




9.3 




<r'' > 


8.7 


11.2 




6.2 



an improved optical potential and more precise data it should be repeated with specific nuclear models. However, at 
this stage the proper choice of a{ro) and the form of the proton (charge) densities are more urgent questions. 

C. Constraints on the optical potential parameters 

The NN amplitudes of eq.([3]) are related to the amplitudes tested in scattering experiments. The latter extrapolated 
to the NN threshold yield scattering lengths and scattering volumes which parametrize the low energy scattering. 
The relation of the experimental lengths and volumes to those required in the optical potential is not direct: first 

- as{E) and ap{E) are strongly energy dependent and one needs these amplitudes for bound particles and second 

- some nuclear many body corrections may arise. An additional difficulty is related to the large number of NN 
partial waves involved, and at this stage of research one can operate only with the spin averaged values. So far, the 
safest method to find the optical potential was to extract the potential parameters from the best fit to the atomic 
data @, In our work we adopt a mixed procedure: the parameters are semi-free, subject to constraints from other 
experiments. In addition, we are guided by the Paris model of NN interactions. 

At the nuclear surface p may interact with quasi-free but bound nucleons. In the Np system the relevant energy is 
negative since both particles are bound and some recoil energy is taken away by the relative motion of the Np pair 
with respect to the residual A — 1 nucleons. Hence, one needs to know the amplitude / below the Np threshold 

/ - asi-Es - Erec) + 3pp'ap{-EB - S„c) (10) 

where Eb is the sum of antiproton and nucleon separation energies and Erec — Prec/'^l^rec is the recoil energy. The 
distribution of the recoil momenta - Prec - is calculable from the Fourier transforms of the antiproton and nucleon 
wave functions \['p(r)(^Ar(r). In practice, such calculations can be done easily in the lightest atoms: deuterium and 
helium [2^. In the "lower" {I = 0) and "upper" (/ = 1) orbitals of ^H p atoms one obtains average values < Erec >= 9 
and 4 MeV respectively. In heavier nuclei such calculations are less reliable but the spectrum of Prec was measured 
with p stopped in Ne streamer chamber [30]. The Prec distribution peaks at « 180 MeV and gives < Erec >~ H 
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FIG. 2: The proton and neutron densities in two Ca isotopes calculated in a HFB shell model have been approximated by the 
two parameter Fermi distribution with constant c and variable o(r). Functions a(r) are plotted to show its asymptotic behavior 
at large distances. The arrows indicate absorption radii for the lower (L), upper (U) and radiochemical (R) experiments, 
correspondingly. 



MeV. This value is used in further calculations. In addition, we estimate a 5 MeV difference of < Erec > in the 
"lower" and the "upper" atomic state. 

The energy dependence of the Np scattering amplitudes in the sub-threshold region may be obtained, to some extent, 
from the ^H, ^H, "^He, "^He atoms since the nucleon separation energies in these nuclei span the region from to 21 
MeV. Calculations based on the multiple scattering series summation from Ref. [1^ and data from Refs. 0, 0j 11] j 
are presented in Fig. ([3]). From the level shifts and widths in these elements one can extract averaged absorptive parts 
of Im 05 and Im ap via a best fit procedure [isj . The nucleon binding energies characteristic for surfaces of large 
nuclei locate the Eb + Erec energies in the sector (40 to 15) MeV. Corresponding values of Im as and Im ap indicate 
the strengths of the absorptive optical potentials expected in nuclei. Two results are of significance in the analysis of 
the optical potential. First, the 5*- wave absorption strength Im as increases with the decreasing energy. The physics 
behind it is related to a broad quasi-bound TV TV state which is indicated by the atomic data and generated by the Paris 
potential. Its existence is also inferred from the pp correlations observed in J/i/; — > ppj decays, [l3j. The impact of 
such phenomenon is discussed in the next section, it is clearly seen in the comparison of "lower" and "upper" widths. 
Second, a resonant-like behavior arises in a P-wave close to -10 MeV. Figure ([3]) indicates that a similar effect is also 
generated by the Paris potential model, where it is attributed to the iso-triplet, spin-singlet ^^Pi quasi-bound state. 
It affects antiproton capture on very loosely bound nucleons p^ . 

A different constraint on the isospin structure of Im 05 follows from studies of n mesons produced by p stopped 
in chambers pTl . [T3 |. These experiments yield ratios of annihilation rates Rn/p = (^{pn) / (t{Pp) given in Table HIl 
Since the p are stopped in light elements the interaction involves S waves and the cross sections are related to the 
absorptive amplitudes a{Np) ~ Imas(iVp) for the np and pp pairs, respectively. Inspection of Fig. ^ shows that the 
Rn/p{S) for the S waves is well reproduced by the Paris potential model. The same model predicts Rn/p{P) ~ 1 for 
the average of P waves at energies of our main interest. These constraints are considered as possible improvements of 
our optical potential. It turns out that the condition i?„/p(P) « 1 is consistent with the atomic data. In this respect 
our results follow the findings of Ref. [l3|. The chamber result Rn/p{S) « 0.5 leads to no or a marginal improvement 



7 



n3 ^ 



n3 ^ 



Im a„ 




-2 



-1 



-40 -30 -20 -10 10 

ENERGY (MeV) 




-1 



^He + np 



-40 -30 -20 -10 10 

ENERGY (MeV) 



FIG. 3: Left panel: the absorptive parts of pN S'-wave amplitudes in the sub-threshold region extracted from the atomic level 
shifts and widths in ^H, ^H, '^He, *He. Im as{E) averaged over spins and isospin of the nucleons are given by squares located 
at the corresponding values of Eb+ < Erec > ■ The curves are calculated with the recent updates of Paris NN potential, [l3|. 
The dotted lines refer to pp and the continuous lines to np systems. Right panel: the absorptive parts of pN P-wave amplitudes 
in the sub-threshold region extracted from the atomic level shifts and widths in '^H, ^H, ^He, ''He (circles) The curves - Paris 
potential calculations. 

over that for i?„/p(S') « 1. On the other hand, it changes the neutron density radii extracted from p atoms. 



TABLE II: The experimental antiproton capture ratios R^/p = a{pn)/a{pp) extracted from capture in flight. 



Element 


-^n / p 


reference 




0.81(3) 




^He 


0.47(4) 


[12j 


'He 


0.48(3) 


[12j 



D. The constrained best-fit potential 

An overall best fit to the 117 atomic data was performed. It includes measurements in N, C, O, P, S, CI, Ca, Fe, 
Co, Ni, Y, Zr, Cd, Mo, Sn, Te, Pb, Bi, Th and U, in total 37 isotopes. Several precise results were excluded: ^'^Li 
- as the optical potential may not be well applicable there, "'^^O - due to the large uncertainty of the charge density 
profile as exemplified already in ^^O, Yb - due to its uncommon deformation. An anomalous lower width in ^^Fe was 
dropped and several very old and uncertain results have not been considered. 

One cannot find the best fit parameters jointly with the details of uncertain neutron densities in each individual 
isotope. To obtain an overall best fit the following strategy was adopted. An overall trend of the neutron rms radius 
as a function of the neutron excess S = {N — Z)/ (N + Z) was assumed and for Z > 10 a linear interpolation was 
tried. Two such trends were fitted before to the atomic p data. One, called "hadronic" is indicated by the hadronic 
scattering data [Slj 



ARr. 



Rrmsin) - Rrmsip) - -0.09(2) + 1.46(12) S. 



(11) 



It is close to the trend obtained by Friedman et al as a result of pionic atom data and the best-fit to the antiproton 
X-ray data [l3] ■ Another " atomic" trend follows from the p atomic data analyzed with zero NN annihilation range 



and Rr 



■/p 



1, 



Rrmsin) - Rrmsip) = ^0.03(2) + 0.90(15)5 



(12) 



Both slopes can be supported by nuclear model calculations and it is appropriate at this point to remind that nuclear 
models are not able to predict the neutron radii and we have to fit some parameters to experimental data. This 
point is particularly strongly stressed by Furnstahl (ssj and a similar point of view is taken in the recent extensive 
calculation by Ring et al. [27\ . Neither of these trends indicated above needs to be true, some nuclei lie far away from 
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the averages given above. Here, we use the slope to get an initial insight into the best fit possibilities. Next with 
the best parameters for the optical potential, we find the best results for the neutron excess in some nuclei. This is 
done, following the procedure of Ref. [l^], in terms of free neutron diffuseness a„ (neutron halo, c„ — Cp) or neutron 
half-density parameter c„ (neutron skin, a„ — ap). An additional question is whether the inclusion of extra data, the 
radiochemical measurements or other experiments, determines some correlation of these two parameters. 

TABLE III: Overall parameters for the optical potential fitted to atomic data: 78 - lower level shifts and widths, 39 - upper 
level widths. In the first line the slope follows eq |12l With other "hadronic" cases the slope was varied but close to the one 
given by eq llll The best result is obtained with AR'^p — —0.10(2) + 1.65(5)5. The r^ms denotes the root mean square radius 
of the folding Gaussian profile. All charge densities come from Ref. _24]. Factor fi =4/3 for protons and 2/3 for neutrons 
describes the chamber result, a„ or Cn denotes the neutron density free parameter. 



"slope" 


parameter 


data 






as 


[fm] 


ap 


[fm«] 


rrms [fm] 


Eq[l2] 


an 


aU 


293 


2.50 


-1.70 


- 0.95i 


-.01 


- 0.391 


0.80 


AR'^p 


Cn 


aU 


285 


2.43 


-1.70 


- 0.951 


-0.0 


- 0.351 


0.79 


ARnp 


an 


aU 


262 


2.24 


-1.65 


- 0.85i 


-0.0 


- 0.411 


0.79 


ARnp 


an 


aU 


262 


2.24 


-1.60 


-/,0.95i 


-0.0 


- 0.391 


0.79 


AR„p 


an 


lower 


160 


2.24 


-1.60 


- 0.951 


-0.0 


- 0.411 


0.79 


ARnp 


an 


upper 


87 


2.05 


-1.60 


- 0.621 


-0.0 


- 0.411 


0.79 



The initial best fit results are given in Table IIIII The faster ARrms slope deduced from the hadronic scattering 
offers a better x^- This result differs from the result obtained in ref., [1^, as it is based on both widths and level 
shifts while ref. , [11], used level widths. The best result Ai?„p = —0.10(2) + 1.65(5)(S obtained with a change of 
the diffuseness parameter a„ corresponds to line 4 in Table III. The difference of our result and that given by eq. (fTTj) 
and Ref. [lOj] is due mainly to the input charge densities. Here, the more recent muon-based data |24] are used while 
results of Ref. are based mainly on the electron data. As discussed in previous sections and indicated in Table I, 
the muonic data are close to the "average" results and in this sense seem to be preferable. In addition. 

Another minimum of is obtained with changes of the half-density radius c„ at a very high slope ARj-ms — 
—0.10 -I- 2.0 S. Although the total is much worse in this case, it turns out to be better for some specific isotopes. 
This question is discussed in the next section. Let us also notice that the best potential parameters are fairly close 
for both the "atomic" and the "hadronic" slopes. 

The preference found in chamber experiments [Tll.[l^ for a stronger S wave annihilation on protons, i?„/p(S') « 0.5, 
is now introduced into the optical potential parameters. The effect is given in line 4 of Table Hill It does not change 
the best fit, but is accepted by the data. It is also weakly reflected in the analysis of radiochemical data as the 
absorption in P waves is dominant. The effective ratio of the absorption rates Rn/p depends on the partial wave in 
the antiproton-nucleon system. For each atomic state the mixture of S and P wave is slightly different. We calculate 
a kind of average value characteristic for all Z, n, L states considered in this analysis. For the best fit potential of 
table [TVl one obtains Rn/p = 0.86(4) as the best average representation for the capture ratio a{pn)/a{pp) at large 
distances tested in the radiochemical experiments. The uncertainty given in parentheses describes the dispersion of 
Rn/p obtained in this way. A more detailed discussion of this point may be found in Ref. [34] . 

Lines 5 and 6 in Table [1111 give separate best fits to the lower and upper levels. To find the best solution on those 
limited data sets only Im as was varied. The motivation for such a choice comes from the shape Im as{E) in Fig. ([3]). 
The best fit absorptive parts compare well with the values of Im as{E) at the characteristic subthreshold energies of 
about -30 MeV. 

TABLE IV: The best fit potential based on X-ray data, consistent with the chamber experiments, the lightest atoms and 
the N — N Paris potential, fi = 4/3 for protons and 2/3 for neutrons. as{"upper") Is to be used for the upper levels and 
as(" lower") for the lower ones. 



x' 


xVn 


as {" upper" )[fm] 


as {" lower" )[fm] 


ap [fm^] 


ro[fm] 


247 


2.11 


-1.60 -1 0.74/i 


-1.60 - 1 1.10/i 


-0.0 - 10.39 


0.79 



Certain incompatibility in the description of lower and upper widths was noticed already in Ref. [9] . This effect is 
also reproduced here and it is indicated in Table IIIII The upper levels require weaker absorption and the explanation 
comes from the left panel of Fig. ^ . The p in the upper level encounters less-bound nucleons and the recoil energy is 
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also smaller. As discussed above, one expects the difference in recoil of some 5 MeV and the difference in the average 
bindings of some 3 MeV. The central energy involved in as,p{—EB — E^ec) amounts to about - 35 MeV in the lower 
levels and about - 27 MeV in the upper levels. The total 8 MeV energy shift in the argument of Im as{E) may reduce 
its value by 0.4 fm as required by the X-ray data. Table Hill indicates that such a change reduces the total from 262 
to 247. The final results are summarized in Table [TVl This potential does not offer the best fit to the data, however 
it offers the best fit under the chamber constraint on the S wave absorption Rn/p{S) = 0.5. Relaxing this condition 
one can improve the total by 2-3 units. In addition, it is possible to improve the fit to the data by relaxing the 
other condition i?„/p(P) = 1. There is no direct experimental indication for this P wave constraint, we are motivated 
entirely by the Paris potential calculations. 

III. THE NEUTRON RADII 

The differences of neutron and proton mean square radii, extracted from from several atoms, are given in Tables fVl 
and IVII These results indicate two basic problems: 

First, the extracted Rrms{n) — Rrmsijp) depend rather strongly on the charge density input. In particular, the 
results in the Sn isotopes depend on the < r^^ > moments dominating the upper level width. For two charge profiles 
jj^ 112, 116,120, i24gj-^ . Qj-^g from Ref. [2^1 (electron scattering) and the other from Ref. ^ atoms) the ratios < r^^ >e 
/ < r^^ >^ are : 1.63, 1.32, 1.42, 1.09, respectively. These differences reflect on the differences in neutron radii 
extracted from the antiproton data. As discussed above, the charge densities given by Fricke et al [2^ are close to 
"average" densities generated in several /i atom and electron scattering experiments. In this sense these results are 
more likely than the others. 

Second, the minimum x^ in each isotope was obtained either by enhancing the diffuseness parameter a„ or half 
density radius c„. In most cases both minima offer good Xpd/ ~ 1 additional data are required to determine the 
nuclear surface shape. However, in ^^^Sn, ^^^Sn, ^"Zr and ^"^Pb the fit is bad and Xpd/ ~ ^- ^^'^ discrepancy comes 
from the level shifts. In these cases the change of the half-density radius offers better x^ and such case is given in 
Table IVII for Pb. The neutron radii obtained via the c„ extension are excessively large and run into conflict with 
most of the other data. We exemplify this situation in the case of Sn and Pb atoms, but it is typical to other large 
nuclei. In the Pb nucleus the two minima for Rrms{n) — Rrms{p) yield very different results. The solution obtained 
be enlarging a„ is close to the results obtained in (p,p') scattering experiments Presumably it is the one that is 
physically acceptable. The uncertainty related to the charge density input is rather large but again there are good 
reasons to favor the last column based on the Ref. [24| . It is interesting to note that the second solution characterized 
by the extension of c„ and large difference Rrms{n) — Rrms(p)of 0.5 fm is close to the result first obtained in the 
neutron/proton pickup experiment of Koerner and Schiffer [36|. This solution, if it represents the reality, has an 
interesting astrophysical significance [33| . 

The same interpretation follows from the results obtained in Sn. The neutron radii obtained with enhanced a„ are 
close to those obtained in the proton scattering experiment [13] that yiel ds Rr,ns{n) — Rrms{p) — 0.25(5) fm in ^^'^Sn, 
larger than 0.18(7) extracted via the dipole-state excitation method [S^. The enhancement of c„, allowed by the x^j 
yields very large neutron radii not conflrmed in other experiments. 

It was argued in Ref. [2^ that the additional data needed to pinpoint the nuclear surface shape come from the 
radiochemical measurements. With an optical potential derived from zero range NN interactions, these data favored 
the halo type solution. Such conclusion is also supported here, and exemplifled in Table IVl by the ^^''Sn case. It also 
favors the halo type solution but the extracted Rrms{n) — Rrms{p) are larger and closer to the hadron scattering results. 
The same effect occurs in the case of Pb nucleus. The favored neutron radius excess Rrms{n) — Rrms(p) = 0.22(3)fm 
given in Table IVTl is 0.05 fm larger than the radius obtained with different optical potentials in Ref. [331 ■ There are 
two factors contributing to that: the first and the dominant one is the chamber data that enforces smaller i?„/p, and 
second the a(r) given by nuclear models tends to be smaller at large distances. Due to both factors the neutron radius 
obtained from the X-ray data is slightly smaller than the radius obtained from the joint X-ray and radiochemical 
data. 



IV. CONCLUSIONS 

In this paper a step was taken to describe the antiprotonic atom data in a semi-phenomenological way. The best-fit 
optical potential was constrained by the atomic data from the p hydrogen, deuterium and helium atoms. Additional 
information on the S wave isospin structure was extracted from the chamber low-energy p data. All these constraints 
refer to the absorptive potentials. Following results have been obtained: 
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TABLE V: The Rrms(n) — Rrms{p), in fm units, extracted from the X-ray data in antiprotonic Sn, Te, Zr and Ca atoms. The 
first column refers to charge density profile the second indicates the free parameter of the neutron density. As discussed in the 
text, the results obtained with charge density from Ref. [SJ are the favored ones. The last column and the entries marked by 
* are extracted from the X-ray and radiochemical data. 







"'Sn "'^Sn ^20g^ 124g^ 


i24sn* 


[24.1 

2pF[23] 
2pF[23| 


In 

a„ 

C,i 


0.21(5) 0.22(4) 0.22(4) 0.23(4) 
0.12(4) 0.15(4) 0.16(4) 0.21(4) 
0.25(5) 0.31(6) 0.35(8) 0.42(7) 


0.26(4) 
0.25(3) 
0.45(6) 






122^^ 124^g 126^6 128Te 


128rpg* 


[40] 


ffln 


0.10(7) 0.04(4) 0.12(5) 0.08^;°,! 


0.15(5) 






"«Te i30Te* s^Zr ''''Zr 


s«Zr* 


[40], [24] 


ffln 


0.10(5) 0.17(4) 0.08l;!34 0.12(3) 


16+006 

'-'•^^"-0.08 






«'Ca 42 44 c^ 48 


«Ca* 


[24] 


ffln 


-0.09(9) 0.01(6) .02(6) 0.09t;o8 


0.10(7) 



TABLE VI: The 7?rms(n-) — Rrms{p) differences extracted from antiprotonic Pb atoms. The entries in parentheses [ ] 
give the corresponding given mostly by the level shifts. Second line gives the difference of neutron and proton diffuseness 
parameters a„p = a„ — ap 



parameter 


Pcharge l'23] 


Pcharge f^^j 


On 


0.16(3) [10] 


0.22(3) [10] 


dnp 


0.11 


0.15 


Cn 


0.47(8) [3] 


0.55(8) [3] 



• A consistent description of the annihilation parameters in antiprotonic X-ray data, chamber data and the 
updated Paris potential model, which incorporates the NN scattering. 

• This consistency allows to separate the S and P-wave contributions to the optical potential for antiprotons and 
estimate fairly precisely the ratio of and pp annihilation at very distant nuclear surfaces. The corresponding 
parameter i?„/p = 0.86(4) averaged over a range of nuclei is obtained at these distances and it allows to discuss 
jointly the X-ray data and the radiochemical data. 

• A definite energy dependence of the NN scattering amplitudes in the subthreshold region is indicated by the 
lightest atom data, chamber data and the NN Paris potential. For the S wave it is reflected in the widths of 
upper and lower levels in heavy atoms. For P waves there is an indication of a fairly narrow NN quasi-bound 
state. It is likely that such a state has sizable effect in a small sector of the radiochemical data taken on loosely 
bound protons. On the other hand, apart from the deuteron, the X-ray data offer no convincing evidence for 
such a state. New X-ray experiments performed with nuclei of small nucleon separation energies would be 
helpful to resolve this question. 

• The p X-ray data, although related to the nuclear surface, may supply important information on the neutron 
radii. However, these data taken by itself cannot tell precisely if the neutron excess forms an extended half- 
density radius, enlarged diffuseness or some specific correlation of both parameters. Jointly with the scattering 
and/or radiochemical data the p results favor the neutron profiles of enlarged diffuseness. 

• One important source of uncertainty in the description of the p nucleus interactions is the real part of the optical 
potential and the related level shifts. It is complicated by the increasing evidence of pN quasi-bound states. 
Precise measurements of the atomic fine structure, in particular in very light atoms, would be very helpful to 
resolve this question. 
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APPENDIX A: THE GRADIENT POTENTIAL 



The P wave colhsions of N and N, described in their center of mass system, involve the relative momentum 
P = (Pn — Pn)/2- It generates the gradient potential given by expression Q in the main text which contains the 
nucleon wave functions and their gradients. These are not easily reducible to the nuclear densities. Summation over 
nucleon states and calculations of the derivatives in eq.Q yield involved expressions for Vq- The input uncertainties 
and experimental errors call for a simpler result. To obtain it let us remind that in the "upper" high angular 
momentum states \n,l > the level shifts e and widths T are quite accurately given by 

e-ir/2~< n,Z|F°P*|n,/ > . (Al) 

For such average values the gradient potential Vq may be expressed by a simpler formula '4l'| 

<nJ\VG\n,l>^^^ap f dngp{u) f drp(r - u)^^*^ ,(r)+ | vI/„^Kr - u) p E^Z^V'W + (A2) 

where p = E | (/Sq, p is the nuclear density, '^n.i is the atomic wave function and 

DHl,{r) = iK'nAr) + iii±llvl;2^^^^(r-)]i-, (A3) 



D'^lir) = l^U-) + ^^^^^LMlf • (A4) 

In these equations '^r,n,i denotes the radial part of the atomic wave functions in states of main quantum numbers n 
and angular momentum I. The nucleon angular momentum is denoted by L. Eqs. (jA3p and (|A4p are obtained with the 
"gradient formula" which splits the derivative into tangential and radial components. The two terms in eq.([5]) which 
contain mixed n gradients generate y™'^ . The latter may be expressed by a formula obtained in Ref. [4l[ , 

< Vm^x / / ^yj^^*;^^^^^^^^ ^ + TrK"^ x)g.yPy + i'^H^ N N?Tx9xyTv] (A5) 

P^NN ^ J J II 

A shorthand notation is used and arguments were put into indices i.e. Qxy — ^(x — y), ^ = ^(x) etc. The first and 
the second term in this equation tend to cancel strongly. The last term, where j denotes the nucleon and antiproton 
currents yields the main effect. For even A (spin zero) nuclei the tangential currents average to zero and one is left 
with the radial components which generate vE'r,™,;^'^ „ iVr,a^r,a contributions. However, in odd-A nuclei the tangential 
components in the last term introduce correlations of the atomic and nuclear currents. It leads to splitting of the 
atomic levels. Usually, this effect is small, but it may be magnified if both antiproton and the odd valence nucleon 
have high angular momenta. 

Now some approximations are introduced to make the gradient potential applicable to practical calculations. For 
nucleons, the radial gradients are needed only at large distances. In this region one has, on average, iy9,. ~ Lpr/'i.a where 
a is the surface thickness parameter. To account for the tangential gradient, the angular momentum factor L{L + 1) 
in ea. (jA3[) is averaged over the three uppermost nucleon shells given by simple shell model [i^. The gradients of 
atomic wave functions may be given explicitly for the circular orbits of interest. These functions are 

*r,«,i(0 = iV(n,Z)r'exp(-r/B7i)F„„,,(r), (A6) 

where N is a normalization, B is the Bohr radius and Fn^ci describes the deformation of the Coulomb wave function 
due to short range (nuclear + finite charge) interactions. The dominant effect in F conies from the damping due to 
absorption. For the radial derivative one has 5'j,„ ; = (V^ ^ 1/nB + F' / F)'^r,n,i and the last term is calculated in 
a quasi-classical way in terms of pwkb{V) — ■\/ 2A/^ [E — V) the value of local momentum of the antiproton inside 
the nucleus. One has F' /F = ipw k b {ycentr + Vcui) — ipwRBiVcentr + Vcui + T^°''*) and numerically one finds that 
dropping the F' / F term altogether yields a minute change of the best fit parameters, making an overall worse by 
2.5. Also, it was found that a numerical calculation of d'^ /dr makes no substantial improvement. 

The approximation (|A2p is well fulfilled in the "upper" atomic orbits. Thus equation (|A4p is used here to define 
the local equivalent to the nonlocal gradient Vg potential that gives the same expectation value in a given atomic 
orbital of a given atom. Expressions (jA2IA3IA4[) lead to a typical folded optical potential 

Vp{r) = ^ap^ / dnfp{xi)Dlp{v - u) (A7) 

p.NN 4 J 
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where the effect of antiproton and nucleon momenta is now included into a function 

r nB 4a^ r nB a 

These approximations bring the optical potential back to the form given by the basic equation [T] in the main text. 
Now Vp is given by eg I ATI summed over protons and neutrons. The interpretation of three terms in ea lASI is fairly 
transparent. The first piece contains tangential and radial momenta of the orbital antiproton. It is the "localized" 
version of the gradient term used in mesonic atoms. The second term contains radial and tangential components 
of the nucleon momenta. The radial ones are expressed in terms of the nuclear density diffuseness parameter. The 
last piece is the mixed term which contains only the radial term as for spin zero systems the average product of 
tangential momenta vanishes. All together the nucleon momentum part contributes about one quarter of On 
a phenomenological level it may be approximately included into the Vs potential term. The difference arises when 
one attempts to relate as and ap to the scattering data. Another special effect of this nucleon term is due to its 
dependence on the nucleon angular momenta. Some enhancement arises in high L shells e.g. in Pb nucleus (12/i 
protons and 14z neutrons in the valence shell). Unfortunately, the L(L+ 1) averaging procedure contains some model 
dependence but the effect is moderate anyway, since the < L{L + 1) > term constitutes less than 10% of the total 
D^. In practical calculations we used a smooth, approximate interpolation < L{L + 1) >= 2 + (3Z/2 — 10)/4. Rather 
slow but systematic increase of follows the rising atomic number Z . It is due to the increase in the L and I values. 
The impact of the nucleon < L{L + 1) > term is also rather limited. Dropping it requires a change of ap from the 
value 0.0 - ^0.41/771^ (line 3 in Table III) to the value 0.0 - iQAif-nr' and the increases by 4 units. 
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